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Abstract 

We describe a measurement of the time-integrated luminosity of the data collected by the BaBar experiment at 
the PEP-II asymmetric-energy e + e~ collider at the T(4S), Y(3S), and T(2S) resonances and in a continuum region 
below each resonance. We measure the time-integrated luminosity by counting e + e~ — > e + e~ and (for the T(4S) only) 
e + e~ — > yfyT candidate events, allowing additional photons in the final state. We use data-corrected simulation to 
determine the cross sections and reconstruction efficiencies for these processes, as well as the major backgrounds. 
Due to the large cross sections of e + e~ — > e + e~ and e + e~ — > yfpT , the statistical uncertainties of the measurement 
are substantially smaller than the systematic uncertainties. The dominant systematic uncertainties are due to observed 
differences between data and simulation, as well as uncertainties on the cross sections. For data collected on the T(3S) 
and T(2S) resonances, an additional uncertainty arises due to T -> e + e~ background. For data collected off the T 
resonances, we estimate an additional uncertainty due to the time-dependent efficiency variations, which can affect 
the short off-resonance runs. 

Keywords: BaBar integrated luminosity 
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1. Introduction 

The BaBar detector operated at the PEP-II 
asymmetric-energy <? + e~ collider and collected physics 
data from October 1999 until March 2008. Most of the 
data were collected at an e + e~ center-of-mass (CM) 
energy corresponding to the mass of the T(4S) 
resonance 01 ■ This "on-resonance" Y(4S) sample 
contains (464.8 + 2.8) x 10 6 BB events $M and is 
used for the study of B-meson decays, CP violation, 
and B° - B° mixing. Data samples collected at the 
T(3S) and Y(2S) resonances in 2008 are used for 
bottomonium studies and for dedicated new-physics 
searches. For each Y(nS) resonance (n = 2,3,4), an 
"off-resonance" sample was collected for studying 
continuum e + e~ — > qq events, where q is a u, d, s, 
or c quark. The off-T(45) sample has a CM energy 
about 40 MeV below the Y(4S) peak mass, and the 
off-T(3S) and off-T(2S) samples are 30 MeV below 
the respective peaks. All on- and off-resonance samples 
are used for charm, t, two-photon, and QCD physics 
analyses. 

Measurements of production cross sections and 
branching fractions often depend on knowledge of the 
time-integrated luminosity X of the collected data sam- 
ple. In some cases, the uncertainty on X is one of the 
major sources of systematic uncertainty |3J. In addi- 
tion, in T-resonance data analyses, background char- 
acteristics or the level of continuum background con- 
tamination are often determined from the off-resonance 
sample. This requires knowledge of the ratio of the 
integrated luminosities of the on-resonance and off- 
resonance samples. 

In this article, we describe the final analysis of the in- 
tegrated luminosity of the dataset collected by BaBar, 
incorporating the latest processing and reconstruction 
of the dataset, improved techniques, and reduced sys- 
tematic uncertainties relative to previous measurements. 
The integrated luminosity is measured with Bhabha 
(<? + <?~ — > e + e~) and dimuon (<? + <?~ — > ju + /lT) events. 
These processes have large, well-known cross sections 
and simple signatures that are easily identified, thus en- 
suring high signal-to-background ratios. We use dipho- 
ton (e + e~ — > yy) events to estimate some systematic un- 
certainties and in the determination of the Y(2S, 35) — > 
e + e~ background contamination. We do not use dipho- 
ton events to directly measure the integrated luminosity, 
due to the significant uncertainty on the cross section 
for this process, as calculated by available Monte Carlo 
(MC) generators. 

The analysis technique and results are presented here 
as a resource for future BaBar physics publications, as 



well as future integrated-luminosity measurements at 
other e + e~ colliders. 

2. Detector and Dataset 

The BaBar detector is described in detail in Ref. [01, 
and only a brief description is given here. Charged- 
particle trajectories are measured with a five-layer sil- 
icon vertex tracker and a 40-layer drift chamber (DCH) 
in a nearly uniform 1.5 T magnetic field. Charged 
hadron identification is provided by a Cherenkov detec- 
tor, and photons and electrons are detected in a CsI(Tl) 
electromagnetic calorimeter (EMC). Muons are identi- 
fied with resistive plate chambers and limited streamer 
tubes inserted between the iron layers of the magnetic- 
field instrumented flux return (IFR). 

A two-level trigger system, composed of a hardware 
("level- 1") stage and a subsequent software ("level-3") 
stage, is used to decide whether an event is recorded. 
Both trigger levels use information from the DCH 
and EMC and employ fast EMC-cluster and track- 
reconstruction algorithms. IFR information is also used 
in level 1. Events passing the level- 1 and level-3 trig- 
ger selections are recorded. After additional prescaling 
(discussed below), events are processed by the offline 
reconstruction, where more sophisticated algorithms 
use information from all detector subsystems. After ini- 
tial stages of the offline reconstruction, an event selec- 
tion and classification stage referred to as the offline fil- 
ter takes place. Classifications of the level-3 trigger and 
the offline filter are used to preselect events for subse- 
quent data analysis. 

The integrated luminosity and its uncertainties are 
determined separately for several data samples. The 
Y(4S) sample is divided into six runs, labeled Run 1 
through Run 6. Each run corresponds to a data-taking 
period with typical shutdowns of no more than a few 
days or weeks. Shutdown periods between runs are typ- 
ically several months long. For each run there is also 
an off-resonance sample, collected during short periods 
interleaved with on-resonance data-taking periods. The 
Run-7 sample contains the Y{3S) and Y(2S) data, as 
well as the corresponding off-resonance samples. Run 7 
also includes a dataset collected at CM energies above 
the Y(4S) resonance, which is not included in this anal- 
ysis. Table [TJ lists the data-taking period and Y reso- 
nance for each run. 

To calculate cross sections and detector efficien- 
cies, we make use of simulated MC samples. The 
BHWIDE |H] MC generator is used to simulate Bhabha 
events, and the KKMC [6] generator with the modifica- 
tions described in Ref. U7[] is used for dimuon events. 
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Table 1 : Data-taking period and the resonance corresponding to the 
PEP-II CM energy V* f° r eacn °f the BABM runs. 



Resonance 


Run 


Month/Year 


T(4S) 


Run 1 


10/1999 - 10/2000 




Run 2 


02/2001 - 06/2002 




Run 3 


12/2002 - 06/2003 




Run 4 


09/2003 - 07/2004 




Run 5 


04/2005 - 08/2006 




Run 6 


01/2007 - 09/2007 


Y(3S) 


Run 7 


12/2007 - 02/2008 


T(2S) 


Run 7 


02/2008 - 03/2008 



We also use KKMC to study possible background from 
e + e~ — > t + t~ events. The BABAYAGA generator with 
next-to-leading-order corrections is used to estimate 
the Bhabha cross section systematic uncertainty. The 
EvtGen 19J] generator is used for studying the back- 
ground from T(2S ) and Y(3S) decays in Run 7. We use 
the BKQED [10] generator to generate diphoton events. 
Events produced by these MC generators are passed 
through a full detector simulation based on Geant4 ifTHl 
and are reconstructed and analyzed in the same way as 
the data. 

3. Analysis Method 

For Runs 1-6, the integrated luminosity is measured 
with Bhabha (e + e~ — > e + e~) and dimuon (e + e~ — > 
li + [T) events, which may include any number of radi- 
ated photons in the final state. For Run 7, e + e~ — > p + fi~ 
events are not used, due to significant uncertainty as- 
sociated with the contribution of the T — > ytyT back- 
ground. 

For a particular data sample, the integrated luminos- 
ity is measured from 

_ ^Vcand — M)gd ^ 
v vis 

where A^ can d is the number of selected signal candidate 
events, of which Nbgd events are estimated to be back- 
ground. The visible cross section <x vis is given by 

o- yh = j ^e(Q)dQ, (2) 

where dcr/dQ is the theoretical differential cross sec- 
tion and e(Q) the efficiency for reconstructing and se- 
lecting signal events for a given phase-space point Q. 
The methods for obtaining each of these quantities are 
discussed below. 



3.1. Event Selection 

The event-selection criteria are designed to yield 
samples of high-purity Bhabha and dimuon events, with 
two high-momentum charged-particle tracks in the cen- 
tral part of the detector and relatively little energy taken 
up by radiated photons. We have chosen the selec- 
tion criteria so that systematic uncertainties arising from 
data-MC differences of event distributions are kept to 
a minimum. Electron vs. muon identification relies 
on comparison of the track momentum with the cor- 
responding energy deposited in the EMC. Event selec- 
tion is performed in two steps: preselection, which takes 
place at the level-3 trigger and during offline reconstruc- 
tion, and is described in Section 13.1.11 and final event 
selection, which is described in Section l3.1.2l 

As a basic requirement for tracks at both selection 
steps, the point of closest approach of the track to the in- 
coming PEP-II beams is required to be less than 1 .5 cm 
in the radial direction (r) and less than 10 cm in the 
beam direction (z). 

3.1.1. Preselection 

Tracks used for the level-3 Bhabha event selection 
must have laboratory-frame polar-angle values between 
0.9 rad and 2.5 rad. Most Bhabha events are selected 
by finding two oppositely charged tracks with CM mo- 
menta above 2.0 GeV/c, where at least one of the tracks 
is associated with an EMC cluster with CM energy of 
at least 2.5 GeV. The CM momenta, polar angles, and 
azimuthal angles of the two tracks are required to sat- 
isfy pi + p 2 > 7GeV/c, |f?i + 62 - n\ < 0.5 rad, and 
l^i - <pi — tt| < 0.3 rad. To maintain high efficiency, 
the level-3 Bhabha selection also accepts events with 
a single track, provided there is an EMC energy de- 
position in the expected location, opposite the track 
in the CM frame. In this case, the requirements on 
the track momentum, the cluster energy, and the po- 
lar and azimuthal angles of the track and cluster are 
Ptrack + ^duster > 6 GeV, |0track + 6» c i ustel - - tt\ < 0.2 rad, 
and |</>track + ^cluster — Jt\< 0.3 rad, where these quantities 
are evaluated in the CM frame. 

Bhabha events are recorded not only for luminos- 
ity determination, but also for EMC calibration. The 
Bhabha cross section increases steeply with decreasing 
<? + e~ scattering angle. Therefore, a large fraction of 
events in regions of high cross section is discarded in 
order to reduce the rate of events handled by the data- 
acquisition system without significant detrimental im- 
pact on calibration. This is achieved by assigning each 
trigger-selected Bhabha event to one of seven bins ac- 
cording to , the larger of the laboratory-frame polar 
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angles of the two leptons. For each bin i, only one of 
every Nj events is logged, where the "prescale factor" 
Nj increases with flj^* . This results in a sawtooth distri- 
bution of cos 6^* that is nonetheless more uniform than 
the original distribution and more suitable for calibra- 
tion purposes. The prescale factor applied to each saved 
event is later used to recreate the initial |cos 9\ spectrum 
for use in the luminosity determination. 

Dimuon events are passed by the level-3 trigger based 
on a very loose criterion of a single track with transverse 
momentum pj > 0.6 GeV/c (a value further reduced for 
Run 7) or two tracks, each having pj > 0.25 GeV/c. 
This loose selection is possible due to the fact that the 
<? + e~ — > pfyT cross section is much lower than the 
<? + e~ — > e + <?~ cross section. At the offline-filter stage, 
dimuon event selection requires two oppositely charged 
tracks. The CM momenta of the higher- momentum and 
lower-momentum tracks must satisfy p\ > 4GeV/c and 
P2 > 2 GeV/c, respectively; the sum of the CM polar 
angles of the tracks is required to satisfy 2.8 < 9\ + 62 < 
3.5 rad; and the sum of the CM energies of the EMC 
clusters associated with the two tracks must be less than 
2GeV. 

The diphoton level-3 trigger selection requires two 
EMC clusters. During Run 1 data collection, the CM 
energy of each cluster was required to be at least 0.35 
of the PEP-II CM energy V^- For Runs 2-7, the re- 
quirement was decreased to 0.3 V*- The sums of the 
polar and azimuthal angles of the clusters must satisfy 
\O\+02-7t\ < ao and \\(f>i -4>2\~A < ao in the CM frame, 
where ao = 0.5 rad for Run 1 and 0. 1 rad for Runs 2- 
7. The trigger is rejected if the event has a track with 
p T > 0.25 GeV/c. 

To facilitate offline checks of simulated trigger ef- 
ficiency, a heavily prescaled, unbiased sample of all 
events satisfying the level- 1 trigger is logged. For cor- 
responding checks of the offline-filter stage, a prescaled 
sample of all logged events is kept regardless of whether 
any offline-filter selection is satisfied. The use of these 
"bypass" samples is discussed in Section [3~3l 

3.1.2. Final Selection 

The Bhabha and dimuon event selections for the 
luminosity analysis impose additional, tighter final- 
selection criteria, relying on event properties obtained 
with the offline reconstruction. 

For Bhabha candidates, the CM polar angles of the 
tracks are required to satisfy |cos 9\ < 0.70 rad for one 
track and |cos6>| < 0.65 rad for the other track. We 
require Pi > 0.75 and P2 > 0.50, where the scaled 
momentum P, = 2/7,/ y7 is twice the ratio of the CM 
momentum p t of track i to the PEP-II CM energy ^[s, 



and the index i = 1 (i = 2) denotes the track with the 
higher (lower) CM momentum. The acolinearity angle 
a, defined as 180° minus the CM angle between the two 
tracks, is required to satisfy a < 30°. We attempt to 
geometrically associate each track with an EMC clus- 
ter and calculate the ratio of the cluster energy to the 
track momentum in the laboratory frame. Denoting the 
higher (lower) ratio with {E/p)n {{E/p)i), we require 
(E/p) H > 0.7 and (E/p) L > 0.4. If only one track is 
associated with a cluster, it must satisfy {E/p) > 0.7. 
Events with no track-cluster association are rejected. 

For dimuon candidates, we require |cos#| < 0.70 rad 
for one track and |cos 0\ < 0.65 rad for the other track, 
Pi > 0.85, P 2 > 0.75, and a < 20°. At least one track 
must have an associated EMC cluster with CM energy 
less than 0.5 GeV. If a cluster is associated to the second 
track, its CM energy is required to be less than 1 GeV. 

Diphoton candidates are selected by requiring events 
with two EMC clusters with energies E\,E<x satisfying 
2EJ > 0.85 and 2E 2 / yfs > 0.75. The CM polar 
angles of the clusters must satisfy |cos 6\ < 0.7 rad for 
one cluster and |cos 6\ < 0.65 rad for the other, and the 
acolinearity angle must be smaller than 10°. If there 
are tracks in the event, the track with the largest CM 
momentum must satisfy Pi < 0.5. 

Hadronic events (e + e~ — > hadrons) are used in the es- 
timation of the T — > e + e~ background. We select such 
events by requiring at least three tracks and a primary 
vertex location consistent with the known beamspot. 
The total energy of tracks and clusters must be greater 

than 0.3 -\fs, and the ratio of the second to the zeroth 

I — 1 

Fox- Wolfram moments [12] is required to be smaller 
than 0.95. The distance between the primary production 
vertex of the tracks in the event and the time-averaged 
beamspot position must be less than 0.5 cm in r and less 
than 6 cm in z. 

Figs. [TJ through |4] show examples of the Bhabha 
and dimuon selection-variable distributions for data and 
simulation. Although in some cases there are visible 
differences between the distributions in data and in sim- 
ulation, the loose selection criteria ensure that these dif- 
ferences have negligible impact on the knowledge of the 
signal efficiency. 

3.2. Background Estimation 

3.2.1. Background Sources Common to All Runs 

The efficiency for e + e~ — > t + t~ — > p* p^v^v^VrVr 
events to pass the dimuon selection is determined us- 
ing MC. We find the fraction of such events in the se- 
lected e + e~ — > [i + p~ candidate sample to be (0.0816 ± 
0.0033)%. The fraction of Bhabha events in the dimuon 
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Figure 1 : Distributions of the scaled CM momentum P, = 2p, / t/s and cosine of the CM polar angle 8, for the higher-momentum (i = 1) and lower- 
momentum = 2) track in candidate e + e~ — > e + e~ events in a fraction of the data (Run 4; solid red histograms) and for simulated e + e~ — » e + e~ 
events (dashed black histograms). The simulation histograms are normalized to the area of the data histograms. The upper two rows of figures 
show the Pj distributions with linear (left) and log (right) vertical scale. In each scaled-momentum plot, the vertical line shows the minimum value 
for events that are retained. When plotting each variable, the selection criteria on all other variables are applied. The |cos#,| (i = 1,2) plots are 
made with |cos 9j\ < 0.7. 
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Figure 2: Distributions of the CM acolinearity angle a, and the higher (lower) laboratory-frame energy-to-momentum ratio E/pu (E/pi) for 
e + e~ — » e + e~ candidates in a fraction of the data (Run 4; solid red histograms) and for simulated e + e~ —> events (dashed black histograms). 
The simulation histograms are normalized to the area of the data histograms. The distributions are shown with linear (left) and log (right) vertical 
scale. In each E/p plot (log-scale a plot), the vertical line shows the minimum (maximum) value for events that are retained. When plotting each 
variable, the selection criteria on all other variables are applied. 
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Figure 3: Distributions of the scaled CM momentum P, = 2p, / t/s and cosine of the CM polar angle 8, for the higher-momentum (i = 1) and lower- 
momentum ((' = 2) track in candidate e + e~ — > events in a fraction of the data (Run 4; solid red histograms) and for simulated e + e~ — » // + //~ 
and — > t + t~ events (dashed black histograms). In the log-scale plots, the dotted blue histograms show the small contribution of e + e~ — » t + t~ 
events to the simulation histograms. The simulation histograms are normalized to the area of the data histograms. The upper two rows of figures 
show the Pj distributions with linear (left) and log (right) vertical scale. In each scaled-momentum plot, the vertical line shows the minimum value 
for events that are retained. When plotting each variable, the selection criteria on all other variables are applied. The |cos0,| ((' = 1,2) plots are 
made with |cos 9i\ < 0.7. 
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Figure 4: Distributions of the CM acolinearity angle a and of the laboratory-frame energies of the higher-energy (£#) and lower-energy (£/,) 
EMC clusters matched to the tracks in candidate e + e~ —> y^iT events in a fraction of the data (Run 4; solid red histograms) and for simulated 
e + e~ — > and e + e~ — > r + r~ events (dashed black histograms). In each log-scale plot, the dotted blue histogram shows the small contribution 
of e + e~ — i t + t~ events to the simulation histogram, and the vertical line shows the maximum value for events that are retained. The simulation 
histograms are normalized to the area of the data histograms. The distributions are shown with linear (left) and log (right) vertical scale. When 
plotting each variable, the selection criteria on all other variables are applied. 
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sample is determined in the same way, and is found to 
be (0.02 + 0.01)%. In both cases, the uncertainties are 
due to MC statistics, and are much larger than those ex- 
pected due to uncertainties on the efficiency or the cross 
sections of the various modes. 

To estimate the background due to cosmic rays or 
beam-gas interactions, we select dimuon candidates 
where the point of closest approach of the tracks to the 
beamline is between 10 cm and 30 cm of the interac- 
tion point in z, and that satisfy all other requirements. 
From this sample, the level of contamination of cos- 
mic events in the dimuon sample is determined to be 
(1.8 ± 0.7) x 10" 5 , which we take to be negligible. 

The background level in the Bhabha sample is much 
smaller than the values listed above for the dimuon sam- 
ple, since the visible cross section for e + <?~ — > <? + e~ is 
an order of magnitude larger than for e + <?~ — > yU + //~. 
Therefore, the background in the Bhabha channel is ne- 
glected. 

3.2.2. Y Background in Run 7 

The on-resonance Run-7 sample contains non- 
negligible contributions from the decays T(2S) — > 
e + e~, Y(3S) — > e + e~ and, to a smaller extent, from cas- 
cade decays such as T(2S) — > n + n~T(\S) — > 7i + n~e + e~ 
or T(2S) yXbj(lP) -> yyT(lS) yye + e -. This 
type of background, which we label as T — » e + e~X, is 
negligible in the T(4S ) samples of Runs 1-6. We deter- 
mine the number of Run-7 T — > e + e~X events from 



(3) 



where Nr is the number of e + e — > T events produced, 
and the visible branching fraction 



Z 



(7 e + e-X) = 
Si(Y e + e-X)e t (r ■ 



e + e~X) 



(4) 



accounts for the branching fraction 2?,-(T — > e + e~X) 
and reconstruction efficiency e,(T — > e + e~X) of each 
process (indicated by the index i) that contributes to 
this background. We obtain S v j s (T — > e + e~X) from 
simulated events, generated with branching fractions 
&i(Y — > e + e~X) based on the measurements compiled 
in the Review of Particle Physics Since S(T(3S) — > 
e + e~) has not been measured, we take its value to be 
identical to S(T(3S) — > relying on lepton uni- 

versality in electromagnetic interactions. 

To determine the number of T mesons produced in 
the Y(2S) or T(3S) on-resonance sample, we count the 
number of on-resonance hadronic events and subtract 
the number of off-resonance events scaled by the ratios 



of luminosities and cross sections between the on- and 
off-resonance samples. The luminosity ratio is deter- 
mined from diphoton events. The number of T mesons 
is |0] 



Nr 



off iy yy 



IVyy ) 



1 

Chad 



(5) 



where A^had (N^ d ) is the number of events satisfying the 
e + e~ — > hadrons selection criteria in the on-resonance 
(off-resonance) sample, N y7 (N°y) is the number of 
events satisfying the e + e~ — > yy selection criteria in 
the on-resonance (off-resonance) sample, ehad is the re- 
construction efficiency for the on-resonance hadronic 
events, and a: is a correction factor accounting for the 
small .s-dependence of the visible cross sections of the 
continuum hadronic and yy events. 

Using Eq. (0, we determine that T — > e + e~ back- 
ground constitutes (1.4 ± 0.1)% of the events passing 
the e + e~ — > e + e~ selection in the on-resonance T(2S) 
sample and (0.9 + 0.1)% in the Y(3S) sample. The 
uncertainties are dominated by the uncertainties on the 
T — > e + e~ branching fractions. The uncertainty on Nr 
is 0.9%, dominated by the determination of eh a d, and has 
a negligible effect on the Nr^ e +e- uncertainty. 

In the dimuon channel, T — > yfyT events consti- 
tute (21.9 + 2.2)% of the selected e + e~ — > y. + [T candi- 
date events for the T{2S) sample and (14.3 + 1.4)% for 
the T(3S) sample. Due to the large uncertainty intro- 
duced by this background, dimuon events are not used 
for Run 7, as mentioned above. 

3.3. Visible Cross Sections 

The visible cross sections cr v i s (see Eq. (f2]i) for 
Bhabha and dimuon events are initially obtained from 
the MC simulation for each run period and CM energy^. 
We then correct the values of cr vis for small data-MC ef- 
ficiency differences, determined as follows. 

We determine the inefficiency of the trigger and 
offline-filter selection from the fraction of events that 
fail this selection but satisfy the final selection require- 
ments, using event samples that are allowed to bypass 
the level-3 trigger and offline filter. From the ineffi- 
ciency difference between data and MC, we apply run- 
by-run corrections to cr vis of up to 0.3%. 

The track-reconstruction inefficiency is measured 
from the fraction of Bhabha events in which only one 



7 The MC generators are not valid in some parts of phase space, 
in particular for small-angle Bhabha scattering, which is excluded by 
the analysis selection criteria. Therefore, the simulation can be used 
to evaluate the visible cross section, but not the full cross section and 
efficiency separately. 
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track is found. To minimize the non-Bhabha events in 
this sample, one of the tracks must satisfy tight selec- 
tion criteria: 0.95 < P < 1.05, 0.9 < (E/p) < 1.1, and 
|cos 6\ < 0.70 rad. A second track is not found in 0.2% 
of these events. The identification of these one-track 
events as e + e~ — > e + e~ is justified by the observation 
that the highest-energy EMC cluster, other than the clus- 
ter associated with the track, has CM acolinearity with 
respect to the track of no more than about 10° (some 
acolinearity is expected, since the missed track bends 
in the magnetic field), and that the ratio between the 
energy of this cluster to the track momentum peaks at 
1 . From the data-MC inefficiency difference, we apply 
run-dependent corrections to cr vis in the range 0.14%- 
0.27%. 

Table|2]shows the corrected visible cross sections for 
the different PEP-II CM energies. For Runs 1-6, we 
observe a run-to-run variation of ±0.21% (±0.7%) in the 
value of cr v ; s for the Bhabha (dimuon) channel. 

Table 2: Visible cross section <x v j s (see Eq. (2J) f° r tne different data- 
taking periods categorized according to the center-of-mass energy ^fs, 
which was equal to ("On") or just below ("Off") the masses of the Y 
resonances. Results for the T(45) samples are luminosity-averaged 
over Runs 1-6. The uncertainties are systematic and are described in 
SectionfH 



Sample 


crvis (nb) 


e + e — > e + e 


e + e — > 


On Y(4S) 
Off T(45) 
On Y(3S) 
Off T(3S) 
On T(2S) 
Off T(2S) 


6.169 ±0.041 
6.232 ± 0.044 
6.461 ± 0.037 
6.508 ± 0.056 
6.933 ± 0.042 
6.866 ±0.051 


0.4294 ± 0.0023 
0.4333 ± 0.0025 
0.4488 ± 0.0028 
0.4501 ± 0.0040 
0.4802 ± 0.0030 
0.4721 ±0.0036 



4. Systematic Uncertainties 

Table [3] summarizes the systematic uncertainties, 
which are described in detail below. 

For the selection criteria used in this analysis, we find 
that the cross section reported b y B HWIDE is consistent 
with that of the BABAYAGA Qj}] generator to within 
the statistical uncertainty of the comparison, 0.06%. We 
add this uncertainty in quadrature to the BABAYAGA 
theoretical uncertainty of 0.20% lfl3ll to obtain the total 
uncertainty of 0.21%. The uncertainty on the dimuon 
cross sections is taken to be 0.44%, based on Ref. 01 . 

From the data-MC comparisons described in Sec- 
tion !3.3l we estimate an uncertainty of 0.13% (0.20%) 
for the track-reconstruction efficiency for Runs 1-6 



(Run 7), corresponding to approximately half the largest 
correction within these data samples. An uncertainty of 
0.1% is estimated for the trigger and offline-filter effi- 
ciency correction by rounding up the largest of the run- 
dependent statistical uncertainties of this correction. To 
account for differences between the distributions of data 
and MC events in the variables used for event selection, 
we vary the selection requirements over wide ranges 
throughout the tails of the signal-event distributions, and 
repeat the full analysis for each variation. For each run, 
the largest resulting change in £ is taken to be the asso- 
ciated uncertainty, with values ranging between 0.40% 
and 0.60%. 

The luminosity and systematic uncertainties are eval- 
uated for the entire period of data collection for each 
particular run. Use of subsamples within a run may in- 
troduce time-dependent variations in efficiency that are 
not accounted for in the analysis. In particular, off- 
resonance data are collected at relatively rare intervals, 
and could therefore be subject to such time-dependent 
effects. Therefore, we estimate an additional systematic 
uncertainty for the off-resonance luminosity, account- 
ing for tracking -related and EMC-related time varia- 
tion studied using the on-resonance samples. The on- 
resonance data sample for each run is divided into at 
least ten subsamples with luminosities of about 1 to 
2 fbT 1 each. In each subsample i, we calculate the ra- 
tio Xi = £f e IXj 1 ^ of the luminosity values obtained with 
Bhabha and dimuon events. We use the spread in the 
Xi values, after subtraction of the estimated statistical 
component of the spread, to estimate the off-resonance 
luminosity uncertainty associated with the time varia- 
tion of any EMC-related effects. Similarly, we use the 
spread of the ratios e j '£j y of the luminosity values ob- 
tained with Bhabha and diphoton events to estimate the 
uncertainty due to the time variation of tracking-related 
effects. Finally, these two uncertainties are added in 
quadrature. 

The uncertainties on the background subtraction, de- 
scribed in Section 13.21 are propagated to the final un- 
certainty on For Run 7, we estimate an additional 
uncertainty of 0.2% on the signal reconstruction effi- 
ciency, arising from the uncertainty on the laboratory- 
to-CM boost associated with changing the PEP-II en- 
ergy from the T(4S ) to the Y(3S ) and T(2S ). 

Systematic uncertainties from the different sources 
are added in quadrature, separately for each channel 
(e + e~ — > e + e~, e + e~ — > jU + /0> run, and on/off- 
resonance data-taking period. When combining results 
in Section [5] we take into account the following cor- 
relations between systematic uncertainties. The uncer- 
tainties on the track-reconstruction efficiency and on the 
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Table 3: Relative systematic uncertainties on the measured integrated luminosity. 



Source 


Relative uncertainty on _£ (%) 


Theoretical cross section 


0.26 (e + e-), 0.44 (jx+fT) 


Track-reconstruction efficiency 


0.13 (Runs 1-6), 0.20 (Run 7) 


Trigger & offline-filter efficiency 


0.10 


Data-MC differences 


0.40-0.60 


Time dependence 


0. 16-0.46 (Off-resonance) 


Background subtraction 


0.02 (Runs 1-6), 0.10 (Y(3S)), 0.15 (Y(2S)) 


Boost uncertainty 


0.2 (Run 7) 



trigger and offline-filter efficiency are positively corre- 
lated between the two channels. Uncertainties in the 
theoretical cross section, background subtraction, trig- 
ger and offline-filter efficiencies, and selection-criteria 
variation are positively correlated for the different runs, 
as well as for the on-resonance and off-resonance peri- 
ods. 

5. Results 

Table |4] lists the integrated luminosity results for the 
on- and off-resonance samples for each run. The re- 
sults for Runs 1-6 are averaged over the yU + //~ and e + e~ 
channels, accounting for correlated uncertainties. The 
results obtained with the two modes are compatible and 
have similar overall uncertainties, with the yT uncer- 
tainties being somewhat smaller. (As noted in Section|3] 
the Run-7 luminosity is obtained with <? + e~ events only.) 
The ratios between the on-resonance and off-resonance 
integrated luminosities are also given. Table [5] shows a 
run-by-run breakdown of the results for the T(4S ) peri- 
ods. 
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Table 4: The integrated luminosities of the on-resonance (X n) and off-resonance (X ff) data samples recorded at and just below the T resonances, 
and the ratio between the on- and off-resonance integrated luminosities. The first uncertainties are statistical and the second uncertainties are 
systematic. 



Resonance 


Xon(fb^) 


Xoff (fb^ 1 ) 


-£on / -£off 


Y(4S) 


424.18 ± 0.04 ± 1.82 


43.92 ±0.01 ±0.19 


9.658 ± 0.003 ± 0.007 


Y(3S) 


27.96 + 0.03 ±0.16 


2.623 ± 0.008 ± 0.017 


10.66 ± 0.03 ± 0.03 


Y(2S) 


13.60 ±0.02 ±0.09 


1.419 ±0.006 ±0.011 


9.58 ± 0.04 ± 0.04 



Table 5: The on-resonance (Xon) and off-resonance (X ff) integrated luminosities of the individual Y(4S) runs, and the ratio between the on- and 
off-resonance integrated luminosities. The first uncertainties are statistical and the second uncertainties are systematic. 



Run 


XonCfb" 1 ) 


Xoif (AT 1 ) 


-£on / -£off 


1 


20.37 ± 0.01 ± 0.09 


2.564 ± 0.002 ± 0.014 


7.946 ± 0.006 ± 0.027 


2 


61.32 ±0.01 ±0.26 


6.868 ± 0.004 ± 0.034 


8.928 ± 0.006 ± 0.023 


3 


32.28 ±0.01 ±0.13 


2.443 ± 0.003 ± 0.012 


13.213 ±0.015 ±0.037 


4 


99.58 ± 0.02 ± 0.41 


10.016 ± 0.007 ± 0.043 


9.943 ± 0.007 ± 0.012 


5 


132.33 ±0.02 ±0.59 


14.278 ± 0.008 ± 0.066 


9.268 ± 0.005 ± 0.012 


6 


78.31 ±0.02 ±0.35 


7.752 ± 0.006 ± 0.036 


10.102 ±0.008 ±0.013 
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